Signatures of Wigner Localization in Epitaxially Grown Nanowires 



L. H. Kristinsdottir,^ J. C. Cremon,^ H. A. Nilsson,^ H. Q. Xu,^ 
L. Samuelson,^ H. Linke,^ A. Wacker/ and S. M. Reimann^'Q 
(Nanometer Structure Consortium, nmC@LU) 

^Division of Mathematical Physics, Lund University, Box 118, 22100 Lund, Sweden 
^Division of Solid State Physics, Lund University, Box 118, 22100 Lund, Sweden 

(Dated: December 1, 2010) 

It was predicted by Wigner in 1934 that the electron gas will undergo a transition to a crystal- 
lized state when its density is very low. Whereas significant progress has been made towards the 
detection of electronic Wigner states, their clear and direct experimental verification still remains 
a challenge. Here we address signatures of Wigner molecule formation in the transport properties 
of InSb nanowire quantum dot systems, where a few electrons may form localized states depending 
on the size of the dot (i.e. the electron density). By a configuration interaction approach combined 
with an appropriate transport formalism, we are able to predict the transport properties of these 
systems, in excellent agreement with experimental data. We identify specific signatures of Wigner 
state formation, such as the strong suppression of the antiferromagnetic coupling, and are able to 
detect the onset of Wigner localization, both experimentally and theoretically, by studying different 
dot sizes. 
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The transition to a Wigner crystaP can be viewed as 
a contest between the electronic Coulomb repulsion and 
the quantum mechanical kinetic energy. If the Coulomb 
repulsion dominates, the many-particle ground state and 
its excitations resemble a distribution of classical par- 
ticles located in a lattice minimizing the Coulomb en- 
ergy. In the bulk, the transition to a Wigner crystal 
is only expected for extremely dilute systems^, while 
in lower dimensions, or for broken translational invari- 
ance, it becomes accessible at higher densities^ ^. A lot 
of work has focused on finite-sized two-dimensional quan- 
tum dots^ where the crossover from liquid to localized 
states in the transport properties of the nanostructure 
has been addressed^^ For one-dimensional systems, 
localization has been reported in cleaved edge overgrowth 
structureJ^ and for holes in carbon nanotubeJ^. These 
highly correlated one-dimensional systems exhibit a va- 
riety of fascinating features as reviewed recently^^. Here 
we introduce a third system, based on epitaxially grown 
semiconductor nanowires, which allows a straightforward 
application of tunneling spectroscopy compared to the 
rather involved cleaved edge overgrowth structures and 
avoids further complications due to the isospin degree of 
freedom in carbon nanotubes. 

InSb nanowires^^, as used here, allow for the realiza- 
tion of quantum dots, where the electronic confinement 
along the nanowire is established by Schottky barriers to 
gold contact stripes, see Fig. [TJa). Varying the distance 
between the stripes (here: 70 nm and 160 nm) allows for 
the systematic realization of wires with specific length 
and thereby controlled electron densities. For our cal- 
culations we model the nanowire as a hard-wall cylinder 
with the experimental radius 35 nm. The Schottky bar- 
rier at the semiconductor-metal interface creates a stan- 
dard quantum well with a width equal to the contact 
spacing. The Coulomb interaction between the electrons 
is approximated as that in a cylinder embedded in ho- 
mogeneous matter, taking into account the different di- 
elect ric co nstants of the wire and the surrounding ma- 
teria^^. Exact many-particle states in the wire are 
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Figure 1: (a) SEM-image of the InSb nanowire on a Si02 
capped Si substrate, where the quantum dot is defined by 
Schottky barriers of the gold contacts ('source' and 'drain'). 
Calculated electron density in nanowires of lengths 70 nm, 
160 nm, and 300 nm is displayed in panels (b,c,d), respec- 
tively, for the lowest two-electron states (excitation energies 
are given; 'S' stands for singlet and 'T' for triplet). For the 
two-particle ground state the pair-correlated density is shown 
with the position of one electron marked by a black arrow. 



evaluated with the configuration interaction method. 

The results can be understood in terms of two limiting 
cases: a short wire with no electron localization and a 
long wire with Wigner localization^^. 

The first limiting case, where interaction is dominated 
by kinetic energy, can be described by the independent- 
particle shell model. There the two-particle ground state 
is obtained by populating the lowest single-particle level 
with a spin- up and a spin-down electron. Thus the spatial 
electron density follows that of the lowest single-particle 
level and exhibits a peak in the center of the quantum 
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dot. The lowest excited two-particle state is obtained 
by moving one electron to the first excited single-particle 
level at the cost of the level spacing energy Ae. Thus one 
expects the two-particle excitation energy AE2 ~ Ae. 
Furthermore the spin degrees allow for four realizations 
of such an excited two-particle state, which are typically 
split into a triplet and a singlet due to exchange interac- 
tion. 

In the second limiting case, Wigner localization, the 
electrons are localized at different positions along the 
wire, minimizing the Coulomb repulsion. Thus the two- 
particle ground state density exhibits two peaks and a 
minimum in the center of the nanowire segment. As the 
electrons can have arbitrary spin on each site, one has 
four realizations of this configuration, with a minor en- 
ergy split between a singlet and a triplet. Hence, we 
expect a very small AE2 <C As, while further excitations 
are significantly higher in energy and exhibit a different 
spatial distribution of charge. 

At the onset of localization, the electron density is ex- 
pected to resemble two weakly separated peaks in the 
two-particle ground state. The interaction of the elec- 
trons is substantial, without yet dominating the kinetic 
part. Hence the two-particle excitation energy is consid- 
erably lower than the single-particle excitation energy, 
AE2 < As. However, as the two electrons are not yet 
fully crystallized, AE2 is expected to be well above zero. 

Tunneling spectroscopy is a convenient way to study 
ground and excited states in quantum dot systems. Here 
we can use the gold contacts (Fig. [TJa)) as source and 
drain by applying a bias Vsd between both stripes. The 
nanowire is located on a highly doped Si substrate cov- 
ered by an insulating Si02 layer, which allows for appli- 
cation of a back-gate voltage Vbg providing an approxi- 
mately homogeneous shift in energy of all levels in the 
dot. Varying V^d and Vbg provides the characteristic 
charging diagrams (see e.g.^^) displayed in Figs.[2|c) and 
[sjc) at a temperature of 300 mK. Here high differential 
conductance indicates that the electron addition energy 
(affinity) coincides with the chemical potential in either 
of the gates. The diamonds of vanishing conductance 
centered around zero Vsd are the regions of Coulomb 
blockade, where the chemical potentials of both reservoirs 
are above the energy difference between the (A^ — 1)- and 
A'-electron ground state and below the energy difference 
between the N- and (A/' + l)-electron ground state. As no 
further lines of high conductance are found for lower gate 
bias, we assume that the lowest diamond corresponds to 
N = 1. Half the width of this diamond defines the charg- 
ing energy U. 

Based on the calculated many-particle states, elec- 
tron transport is treated within the master equation 
modeP^^^ with tunneling matrix elements calculated as 
m Ref.G^. The results are displayed in Figs. [2|a) and 
[sja) for the respective experimental samples displayed 
in panel (c). We find that all Coulomb diamonds agree 
rather well, which indicates that the radial excitations, 
which are disregarded in our effectively one-dimensional 
model, only become of relevance for higher particle num- 
bers in the dot. 

Now we focus on the excited states and show, that the 
experimental conductance data along with our theoret- 
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Figure 2: Results for an InSb nanowire of length L = 70 nm. 
(a) Simulated differential conductance as a function of bias 
(Vsd) and gate energy Eg. The number of particles in the 
dot, A, is shown in each diamond, (b) A closer look at the 
area marked by a dashed box in panel (a). The conduction 
lines, where tunneling into the N = 1 ground state and first 
excited state sets in, are marked by the symbols ® and d), 
respectively. The corresponding lines for the entering of the 
second electron, where the dot reaches the N — 2 ground 
state and the A = 2 excited state, are marked by d) and 
® symbols, respectively. The separation between these lines 
provides the excitation energies from the A = 1 and N — 2 
ground states, Ae and A^2, respectively, which are depicted 
by arrows, (c) Experimental differential conductance as a 
function of bias (Kd) and gate voltage (Vbg). (d) Experimen- 
tal differential conductance as a function of magnetic field {B) 
and gate voltage (Vbg). 



ical calculations allow for a verification of the Wigner 
localization scenario described above. 

For a 70 nm wire, the 2-electron density along the wire 
is a single peak, see Fig. [ijb). This corresponds to the 
independent-particle shell model as described above. In 
b) we have marked the lines, where the first elec- 



Fig. 2 



tron enters the one-electron ground state and the one- 
electron excited state, by the symbols ® and (D, respec- 
tively. This reflects the level spacing Ae = 12 meV as 
shown by the horizontal arrow. Similarly, starting from 
the one-electron ground state, the second electron enters 
the dot reaching the two-electron ground state and the 
two-electron excited state at lines marked by the d) and 
(4) symbols. The separation between these two lines rep- 
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Figure 3: Results for an InSb nanowire of length L = 160 nm, 
panels as in Fig. |2] 



resents the excitation energy AE2 = 11 meV. The four 
lines, ©-(4), can be observed in the experimental data in 
Fig. [2|c) (this is clearer for negative bias, as the mea- 
surement results in the positive bias region most likely 
suffer from charging of impurity states). From this fig- 
ure, we read A£^2^^ = 15 meV ~ Ae^^P = 16 meV, and 
hence for the sample of length 70 nm, the experimental 
data are in good agreement with the independent-particle 
shell model discussed above. 

Note that there is some discrepancy between theory 
and experiment regarding the value of Ae and AE2 . This 
could be due to bending of energy levels at the interface 
of the wire and the gold contacts (Schottky barriers), 
which makes the wire effectively shorter than the spacing 
of the contacts. Indeed, simulations of a 60 nm wire give 
Ae 16 meV and AE2 15 meV. 

We can quantify the electron-electron interaction 
strength by the energy difference between the two- 
particle ground state and twice the energy of the lowest 
single-particle level (half-width of the N = 1 Coulomb 
diamond). This provides the charging energy U^^^ = 
6.5 meV for the 70 nm sample, as read from Fig. [2|c). 
That is U < Ae, in accordance with the independent- 



particle shell model being valid when kinetic energy dom- 
inates interaction. 

For the 160 nm wire, the 2-electron density in Fig. [ijc) 
resembles two semi- separated peaks, indicating the onset 
of Wigner localization (as also seen in the pair-correlated 
density). In Fig. 3, the lines (D-® can be identified both 
in the simulation and the experiment. The theoretical 
1.0 meV and Ae = 2.8 meV, as in the 
1.0 meV < Ae^^P = 



results give AE2 
experiments we observe A£^2^^ 
3.2meV. Again, this is in agreement with the scenario of 
onset of Wigner localization discussed above. 

Note that if we would neglect the different dielectric 
consant outside the wire, the onset of Wigner localiza- 
tion would first appear at double the actual wire length. 
Hence the screening due to the different dielectric con- 
stants of the wire and the surrounding material is an 
important effect and must be included in the modelling. 

The energy separation between the singlet and the 
triplet two-electron state, the antiferromagnetic coupling, 
can also be manifested by the magnetic field dependence 
of the differential conductance. The Sz = I part of the 
triplet is lowered in energy by a magnetic field with re- 
spect to the singlet state by qiibB^ where /i^ is the Bohr 
magneton. Fig. [sj^d) shows that there is a level cross- 
ing at Across ~ 0.4 T (marked by an arrow). According 
to Ref.^ the electronic ^-factors are around 40 for two 
electrons in the dot. This provides an energy splitting 
^^mag _ gjj^^ Bcross ^ 1 ^leV in full agreement with the 
calculated value for the 160 nm wire. Note that for the 
70 nm wire, the level splitting is no longer linear in the 
high magnetic field, Across ~ 4 T, at which the cross- 
ing appears (marked by an arrow in Fig. [2|d)). Hence 
we cannot apply the same method to find AE^^^ for 
the 70 nm wire, although its result AE^^^ ^10 meV is 
of the correct order of magnitude. The strong suppres- 
sion of this antiferromagnetic coupling between the two 
electrons (by an order of magnitude, while changing the 
length by about a factor of two) is one of the hallmarks 
of the Wigner crystal stat^^^. 

Finally, our theoretical results indicate complete 
Wigner localization for a 300 nm long wire. Fig. [ijd) 
shows that in the two-particle ground state, the electrons 
are strongly localized, i.e. they form a Wigner molecule. 
From Fig. 4[ b) we observe that the conductance line of 
the N = 2 triplet first excited state (d)) has merged 
into the line of the singlet ground state ((?)), as expected: 
There is no difference in the energy of these two states, 
as there should be no difference between the singlet and 
triplet states of two strongly localized particles. More 
precisely we find AE2 = 9.3 /ieV and Ae = 0.84 meV, 
i.e. AE2 <C Ae. Furthermore we find U = 5.7 meV, 
that is Ae <C U. This conforms to Wigner localization 
being present when kinetic energy is strongly dominated 
by interaction. 

Even for the N = 3 ground state the theoretical cal- 
culations suggest the o nset of Wigner localization in a 
300 nm wire, as seen in Fig. 4^c). The small energy dif- 



ference between the three lowest N = 3 states results 
in a broad conduction line, marked by the symbol (D 
Fig. 4[ b). Unfortunately, we could not obtain exper- 



imental data for this length, since for such a long sam- 
ple and low charge densities the effect of disorder is too 
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Figure 4: Simulation of a 300 nm long wire, (a) Charge sta- 
bility diagram, (b) A closer look at the area in the dashed 
box in panel (a). Symbols ®-® as in Fig. |2] The two lowest 
N — 2 states have approximately the same energy, and hence 
the double conduction line of the 160 nm wire (d) and ® in 
Fig. S d) has merged into a single line leading to the N — 2 
Coulomb diamond. The broad conduction line consisting of 
three lines for the three lowest N = 3 states, is marked by 
the symbol (D. (c) Electron density of the four lowest N — 3 
states. 



strong, creating an effective double quantum dot. This 
can be identified in a charge stability diagram as addi- 
tional kinks in the conductance lines that comprise the 
= 1 Coulomb diamond^^. Such kinks are not present 
in the stability diagram for the 160 nm wire shown in 
Fig. [sjc), implying that disorder has no significant effect 
in that case. Also, Coulomb interaction has been shown 
to decrease the effect of Anderson localization^^. How- 
ever the theoretical results demonstrate the prospects of 
our approach, if more efficient gating schemes are devel- 
oped. 

We have demonstrated the transition from the 
independent-particle shell model to Wigner localiza- 
tion with increasing length of a semiconductor nanowire 
sample. While the excitation spectrum follows the 
independent-particle shell model for the 70 nm wire 
{AE2 ~ Ae), the onset of Wigner localization is ob- 
served for the 160 nm wire {AE2 < Ae) and finally our 
simulations show complete Wigner localization in a wire 
of length 300 nm. There the excitation energy of the 
two-particle state is almost negligible and much lower 
than the level spacing, AE2 <C As, and the calculated 
electron density exhibits two peaks. This shows that 
InSb nanowires form a convenient system to investigate 
strongly correlated systems by well established transport 
measurement techniques. 
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Strategic Research (SSF). 
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